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Abstract 

The interaction of charged particles, moving in a uniform magnetic 
field, with a plane polarized gravitational wave is considered using the 
Fokker - Planck - Kolmogorov (FPK) approach. By using a stochas- 
ticity criterion, we determine the exact locations in phase space, where 
resonance overlapping occurs. We investigate the diffusion of orbits 
around each primary resonance of order m by deriving general ana- 
lytical expressions for an effective diffusion coefficient. A solution of 
the corresponding diffusion equation (Fokker-Planck equation) for the 
static case is found. Numerical integration of the full equations of 
motion and subsequent calculation of the diffusion coefficient verifies 
the analytical results. 



1 INTRODUCTION 

Many efforts that have been made to detect gravitational waves gave, so far, 
no convincing evidence that they were actually being seen [1]. This is due 
to the fact that not only their amplitude is very small [2], but it is highly 
possible that some kind of damping mechanism operates on them as they 
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travel through space [3-5]. This damping may originate in the interaction of 
the gravitational wave with the interstellar matter [6,7]. 

In a recent paper [8], which hereafter is referred to as Paper I, the prob- 
lem of the interaction of a charged particle with a gravitational wave, in 
the presence of a uniform magnetic field, has been considered for various 
directions of propagation of the wave with respect to the magnetic field. It 
was found that, in the oblique propagation case the motion of the particle 
becomes chaotic and may be considered as a diffusion in momentum space, 
provided that its initial momentum is sufficiently large. 

In order to address in detail the interaction of charged particles with a 
gravitational wave, one should try to calculate the diffusion rate (in mo- 
mentum space) of the particles which follow chaotic trajectories. This task 
involves the derivation of a Fokker - Planck (FP) type diffusion equation and 
the calculation of the corresponding diffusion coefficient [9]. 

In the present paper we investigate the energy diffusion of charged par- 
ticles in the presence of a uniform magnetic field, B = Boi z , due to their 
non-linear interaction with a linear polarized gravitational wave, propagating 
obliquely with respect to the direction of the magnetic field (20° < 9 < 60°). 
The analysis is carried out in the framework of the weak field theory, consid- 
ering the gravitational wave as a small pertubation in a flat space time. We 
use the Fokker - Planck - Kolmogorov (FPK) approach and refer to the glob- 
ally stochastic regime, where overlapping of many resonances occurs. In a 
partially stochastic regime the FKP approach can not be applied, as the par- 
ticles do not undergo "normal diffusion" (random walk process) but rather 
follow Levy statistics [10]. This statistical approach is possible only after 
deriving general formulas, that hold for every value of the perpendicular en- 
ergy of the charged particle and not just for the highest ones, as it has been 
considered in Paper I. 

The motion of a charged particle in curved spacetime is given, in Hamil- 
tonian formalism [11], by the differential equations 

where 7r M are the generalized momenta (corresponding to the coordinates x^) 
and the "super Hamiltonian" , H, is given by the relation 

H = \<T - eA,) (tt, - eA v ) ee \ (2) 
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(in a system of geometrical units where h — c — G — 1) In Eq.(2) g^ v denotes 
the components of the contravariant metric tensor, which are defined as 

with rf v = diag(l,—l,—l,—l) and \h lxv \ <C 1. A M is the vector potential, 
corresponding to the tensor of the electromagnetic field in a curved spacetime 
F^y. The mass of the particle is taken equal to 1. For the specific form of 
the magnetic field we may take 

A = A 1 = A 3 = , A 2 = Bqx 1 (4) 



2 THE STOCHASTICITY CRITERION 

We consider the case of a charged particle moving in the curved space- 
time background of a linear polarized gravitational wave, which propagates 
obliquely with respect to the direction of a uniform and static, in time, mag- 
netic field, B = Boe z . The non-zero components of the metric tensor are 
presented in Paper I (see references therein) and we normalize lengths and 
time to ^, where Q is the Larmor angular frequency. We, furthermore, elimi- 
nate one degree of freedom from our dynamical system through the canonical 
transformation 

$ 3 = X° - COS0X 3 , 7T 3 = -cos6I 3 

= X°, 7T = J 3 + Jo (5) 

Accordingly, the problem of the motion of a charged particle in a gravitational 
wave is reduced to a two-degrees of freedom dynamical system [8], and the 
"super Hamiltonian" (2), in this case, is written in the form 

1 2 11 + asin 2 9sin{v&) 2 
B — —1?. ; — ^ — 7Ti 

2 3 2 l + asin{vQ) 1 

1 (x 1 ) 2 1 1 + acos 2 6sin(u<d) 

-cos z 9I. 



2 1 - asin(v&) 2 1 + asin(uQ) 
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1 asin29cos6sin(i/0) 
+ 2 1 + asinuQ) * lh (6) 
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where we have set 

= sindx 1 -^ (7) 

In Eq. (6) a is the normalized, dimensionless amplitude of the gravitational 
wave and v = ^ denotes the dimensionless frequency. 

The dynamical system under consideration possesses chaotic regions in 
phase space when a ^ [7,8,12]. In order to examine the transition from 
regular to stochastic motion we use Chirikov's overlap criterion [13,14] to 
obtain the lowest amplitude of the gravitational wave, a t hi, above which the 
dynamical system shows prominent chaotic behaviour. 

We first write the Hamiltonian (6) in action-angle variables through the 
canonical transformation 

x 1 = (2/ 1 ) 1/2 W 1 , TTi = (2I 1 ) 1/2 cos^\ (8) 

and, since a <C 1, we make the approximation 

; ^ N ~ 1 T asin(uQ) (9) 

The resulting Hamiltonian is of the form 

H = H Q + aH lS in(vQ) (10) 

We expand the trigonometric term of the perturbation in a Fourier series [7]. 
After further manipulation, the Hamiltonian (10) is written in the form 

H = -sin 2 6ll-h 



a 

+ 2 



hsin 2 e Jtiyr)sin{t& - v^) 
e=-oo 

oo 

+ Ji(l + cos 2 9) ]T [2J"(vr) + Jeivr^siniitf 1 - z/tf 3 ) 



e=-oo 

oo 



+ cos 2 6sin 2 9I 2 J2 Jt(yr)sin(£-& 1 - 
e=-oo 



3 A 



oo 



+ I 3 sin26cos0 ]T -^(^)sin(W 1 - v^) 



£=—oo 



(11) 



4 



where Je(£) is the Bessel function of order I, r — (2Ii) 1 ^ 2 sin6 is the linear 
momentum along the x-axis and a prime denotes differentiation with respect 
to £ = vr. The perturbation term of the Hamiltonian function, Hi, depends 
on an infinite series of linear combinations of the angles d 1 and $ 3 , a fact 
that leads to resonances. In this case, Chirikov's criterion states that chaos 
appears when the width of a resonance, SIi, becomes larger or equal than 
the distance between two consecutive first order resonances, AIi. 

By a near identity transformation we remove all trigonometric terms from 
Hi, except from the one of order £ — m, which generates the principal 
resonance and corresponds to the family of islands whose width enters in 
the stochasticity criterion [7]. The resulting Hamiltonian contains only the 
integrable part H$ and the dominant term and it is therefore called resonant 
Hamiltonian, H R [15,16]. Performing the canonical transformation 



tfi* = tf 1 ^3 j* = h 

m 



i; = i 3 + -h 

m 



03* = ^3 ^ 

the resonant Hamiltonian is finally written in the form 



(12) 



-sin 2 6(R - -hf - h 
2 v 3 m ' 



a 



hsin 2 OJ m {0 - h(l + cos 2 6)[2J^) + J m (£)] 



v 



cos 2 esin 2 e(i; - — hyJUO 



+ 2(P 3 --Ii)cos 2 e-J m (0 
m v 



sin(md 



l*^ 



(13) 



Since -# 3 is a cyclic coordinate, the corresponding generalized momentum 
^3 will be a constant of the motion, so that the dynamical system becomes 
one-degree of freedom. Hamiltonian (13) describes the motion of a particle 
around each first order resonance. Using the resonant condition 



d 0l d 03 



m- 



v- 



dA dA 

and the fact that H f« \, we find the order m of the dominant resonance 



(14) 



m = v(l + 2Ii) 1/2 sin8 



(15) 
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In this case, m ^ vr and not m ~ ur, which was the case considered in 
Paper I, for I± >> 1. This is because in the present paper we are interested 
in a general formula for the stochasticity threshold, valid for every Ii. The 
distance AI ± , between two consecutive first order resonances is calculated by 
Eq.(13) and the fact that Am = 1, 

while the corresponding resonant width 5I\ is given by [15,16] 

8am 2 



8h = 



[m 2 - v 2 sin 2 9) (I + cos 2 9) J'^ 



v A sin A 9 

+ (4m 2 - v 2 sin 2 9)cos 2 9J m 1 1/2 (17) 



Then, Chirikov's criterion, 8I\ > AJi, reads 

- < 8|(m 2 - v 2 sin 2 9){\ + cos 2 9)J'^ 

+ (4m 2 - v 2 sin 2 9)cos 2 9J m \ (18) 

The above relation is the most general form of the stochasticity criterion and 
holds for any value of I±, v and 9. We see that for 1\ ^> 1 it reduces to 

8is 2 r 2 \4cos 2 9J m + (1 + cos 2 9)J' 1 l n \ > - (19) 

a 

which is the corresponding result of Paper I. In this approximation, we may 
obtain an asymptotic form of the stochasticity threshold, by taking r — > oo 
[8]. We obtain 

« > 0-07 , \... l 2Q (20) 
(vr) 5 / 3 cos 2 9 

We see that the stochasticity threshold is a rapidly decreasing function of 
v and r. Therefore chaotic behaviour will appear, no matter how small 
the amplitude of the gravitational wave might be, provided that the initial 
momentum of the particles is sufficiently large. It is clear that in this case the 
high order principal resonances will overlap. On the other hand, if 1\ is small 
the chaotic behaviour will appear only if the wave amplitude is considerably 
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large and the relative frequency small, leading to the overlapping of the low 
order resonances. 

Following the above argument, in Fig.l we give the stochasticity thresh- 
old, 

otthri as a function of the low order resonance m, for different values 
of the wave propagation angle 9 and v = 1.8. Notice that a t h r < 0.2 and 
decreases rapidly as the order of resonance m increases. 



3 THE FPK APPROACH 

3.1 ANALYTIC RESULTS 

Following the FPK approach [9], a diffusion equation for the energy distribu- 
tion function of particles averaged over the phases, ^(i^t), can be written 
for the system described by the Hamiltonian (13) 

* = mK^wj (21) 

To lowest order in a, Eq. (21) describes a diffusion process in the variable 

2 

I± — at constant J 3 . The actual expression for the diffusion coefficient 
D(Ii) depends on the assumptions for the phase dynamics [9,17]. In the 
random phase approximation, it reduces to the quasilinear result [18,19] 

D(I 1 )=na 2 ^m 2 H 2 *(^) (22) 

which, in our case, reads 

D(h) = na 2 ^ m 2 H\ 5 (m 2 - v 2 sin 2 9 [1 + 2h}) (23) 

m 

Around each principal resonance of order m we may associate an effective 
diffusion coefficient, D m , by averaging D(If) over the region between two 
successive first order resonances 

: rm+l 

D rn =< D{h) >= — - / D(h)dm (24) 

To calculate D m we use the facts that Am = 1 and 

5{m - mp) 

4(/(m))= WW)\ (25) 



7 



where itlq is a simple zero of f(m) which, in this case, is given by Eq.(15) 
[20]. Accordingly, we obtain 



D m =< D(I lm ) >= -na 2 mH 2 m 



(26) 



where l\ m is the value of l\ at each principal resonance of order m, which is 
found, from Eq.(15), to be of the form 



hr 



1 



m 



2 \ v 2 sin 2 6 



(27) 



and H lm corresponds to the perturbation term of the Hamiltonian (13) for 
h = hm- Eq.(26), in terms of I lm , reads 



x 



-%a 2 vsin6(l + 2I lm ) 1/2 

3 n2 
J lm (l + cos 2 6)J m + (- + AI lm )cos 2 9J, l 



(28) 



We use the above relation in order to determine the analytical values of the 
diffusion coefficient, as it holds for any value of the parameters. It is clear that 
the diffusion coefficient scales with the wave amplitude a and, through the 
value of Ii m , with the order of resonance m. The diffusion coefficient reaches 
high values at low order resonances (small Ii m ) when the wave amplitude is 
large. In the opposite case (small a) the diffusion becomes effective in the 
range of high order resonances and, thus, in large I\. In both cases for a 
given a the diffusion increases as the action increases. 



3.2 NUMERICAL RESULTS 

For the sake of numerical simplicity and in order to speed up the numerical 
integrations we investigate the case of low order resonances using a = 0.2 and 
v = 1.8 throughout the whole of our numerical calculations. Since the results 
scale with the amplitude of the wave, a, the calculated diffusion coefficient 
is also expected to describe, at least qualitatively, the diffusive acceleration 
at more realistic values of a. 

In order to verify that diffusion of particles, due to their interaction with 
the gravitational wave, occurs, we follow the orbits of a particle distribution 
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on a surface of section, defined as the surface ux° = 2nir. In Fig. 2 the tt\ 
versus x 1 plot and the time variation of the action I\ of the distribution of 
orbits (N = 1000 with initial I\ = 106.5) are presented for the case 9 = 20°. 

Notice the distortion of the principal resonance (of order m = 9) due to 
the overlapping of the secondary resonances. This is due to the fact that, 
for the parameters used, the wave amplitude is large. Thus the overlapping 
occurs in a small time-scale. The diffusion in energy, in this case, is verified 
from the considerable spread around the initial action value l\ = 106.5. 

The numerical estimation of the diffusion coefficient is based on the in- 
tegration of the Hamilton's equations of motion for a number of particles 
(N = 1000), having the same initial action I\ and uniform angle distribu- 
tion. The local diffusion coefficient is related to the average variations of the 
action I x through the expression [18,19] 

D W * < (A/.) 2 > -2« A/, »» (29) 

where 

< A/, > = £ '»W ; hM (30, 

3=1 iV 

<(A/i)2> = £ m) -mo)? m 

3=1 

We have performed a number of computational runs, varying the initial 
value of h (2 < h < 128) and the propagation angles (20° < 9 < 60°). In 
Fig. 3 the numerical estimated local diffusion coefficients for I\ = 3.0 and 
7.0, corresponding torn = 4 and 6 respectively, and for 9 = 60° are presented. 
Notice that the integration time is short, as for longer times diffusion over 
a large number of harmonics dominates, causing strong variations to the 
estimation of D. The plateau value of D is chosen as the diffusion coefficient 
for the above actions. 

In Fig. 4 the analytical and the numerical effective diffusion coefficient, 
as a function of the action I±, for different angles 9 are shown. Notice that 
the diffusion coefficient depends strongly on the propagation angle. There 
exists a good agreement between the numerically and analytically estimated 
values, indicating a power law dependence of the diffusion coefficient upon 



and 
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the action, of the general form 



D{h) ~ d /* (32) 

where the values (analytically and numerically estimated) of the constant d 
and the index k, with respect to the propagation angle 6 are given in Table 1. 
The relative error between the analytically and numerically estimated values 
of the index k varies from 7% to 24%. 



SOLUTION TO THE DIFFUSION EQUA- 
TION 



We can easily solve the diffusion equation for the static case, i.e. d t T 
Then Eq. (21) becomes 



d/i d/i 



0. 



(33) 



We substitute the diffusion coefficient from Eq. (32) to find the solution 
in a power law form 



Hh) = 



To 

1-k 



(fc-i) 



(34) 



with T constant. We must emphasize that the above solution is valid for 
relatively small values of 1\. In the case of very large energies {1\ ^> 1) 
an analytic solution of the diffusion equation can be found, by considering 
the asymptotic form of the effective diffusion coefficient in the large energies 
approximation. For ^ > 1, we have r — > oo and therefore m ~ ^ = vr. 
Then, the perturbation term of the Hamiltonian (13) reads 



Rl 



h (l + cos 2 6)J^ 



m) 



Acos 9J m (m) 



In this case, the Bessel equation becomes 



J" m {m) 

and the asymptotic expansions 



J'm(m) 



m 



J m (m) ~ 0.45m~ 1/3 , J' m (m) ~ 0.41m" 2/3 



(35) 



(36) 



(37) 
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hold [21]. Therefore, in the large energies approximation, the effective diffu- 
sion coefficient reads 

D m {h) = Al[ /% - Bll /% + Cll 3/6 (38) 

where 

A = 2 1/2 na 2 vsin9 a 2 , B = 2 3/2 Tra 2 vsin8 ab , 

C = 2 l ' 2 a 2 vsin9 b 2 (39) 



and 



_ 0.41 1 l + cos 2 9 1.8 1 cos 2 6 

a ~ 2 5/6 ^5/3 s i n 5/30 ' _ 2V6 i/l/3 ^ ^ 



Accordingly, the diffusion equation, in the static case, reads 

dh 

Air - Blf 1 " + Cl[ 



^ C l -.5/6 R ,9/6 , ^13/6 ( 41 ) 



where c is an integration constant. Evaluation of the integral on the r.h.s of 
Eq. (41) is possible only when Iy 7^ (f ) 3 ^ 2 , for which the energy distribution 
function appears a simple pole of order 2 [22]. In this case, we obtain 




,1/6 
i l 



2/3^ 

7 . r.V6\l 

(42) 



(a - bitn 

7-1/6 , . ,1/6 

+ 3 _ I In * + + 2tan~ 1 1 



ll"-d d 



where d — (| ) 1 / 4 . This result is simplified considerably in the perpendicular 
propagation case, i.e. 9 — |, for which Eq. (41) gives 

^=^A 1/6 (43) 

5 DISCUSSION AND CONCLUSIONS 

We have studied the interaction of a charged particle, with a plane polarized 
gravitational wave propagating obliquely (20° < 9 < 60°) with respect to the 
direction of the ambient uniform magnetic field. 
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On the basis of Hamiltonian pertubation theory, previous work on this 
problem shows that the motion of the particles becomes chaotic [8]. Follow- 
ing this, we have derived analytical expressions for the stochasticity criterion, 
thus determining where, in phase space, resonance overlapping occurs, with- 
out any assumption regarding the values of the action and the propagation 
angle of the wave. 

We have verified that diffusion of the particles in action Ii occurs and 
we have applied the FPK approach, in order to derive analytical general 
expressions for the effective diffusion coefficient. Numerical integration of 
the exact equations of motion for particle distributions with the same initial 
action I\ was also performed for the numerical estimation of the diffusion 
coefficient. 

Both methods (analytical and numerical) revealed a power law depen- 
dence of the diffusion coefficient upon the action Ii giving similar results, 
with small variations, on the power law index. Based on these results a 
steady state solution of the Fokker - Planck diffusion equation was found. 

The diffusion coefficient scales with the wave amplitude a and the order 
of resonance m (and/or through the resonance condition with the action 
For small a the diffusion is effective in high order resonances and thus in suf- 
ficient large actions. Diffusion of particles is present in low order resonances 
(small values of the action) only when the wave amplitude is large. In both 
cases the diffusion is increasing when the action is increasing. 

There is also a strong relation between the diffusion coefficient and the 
propagation angle. As the angle decreases, the diffusion coefficient increases. 
This is due to the fact that the lower the angle, the greater is the amplitude 
of the wave for which stochastic motion occurs, leading to the fact that more 
resonances can overlap. 

In conclusion, we believe that the FPK approach may describe in a good 
approximation the interaction of charged particles with a gravitational wave 
in the framework of the weak field theory, where the gravitational wave is 
just a small pertubation in a flat spacetime. It is clear that more work has to 
be done in the realistic case of a curved spacetime and in the full non-linear 
theory. 
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TABLES 



9 


d 





k 




analytical 


numerical 


analytical 


numerical 


20° 


0.095 


0.037 


2.254 


2.455 


35° 


0.076 


0.091 


2.223 


2.153 


45° 


0.052 


0.025 


2.160 


2.407 


60° 


0.017 


0.020 


2.037 


1.980 



Table Caption 

Table I: The analytically and numerically estimated values of do and the 
index k with respect to the propagation angle 9. 
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FIGURE CAPTIONS 



Figure 1: The stochasticity threshold a t hr versus the order of resonance m 
for different values of the wave propagation angle 9 and v — 1.8 

Figures 2a and 2b: Surface of section plots for the case a = 0.2, v = 1.8, 
9 = 20°. N = 1000 orbits with initial h = 106.5 are presented: (a) The m 
versus x 1 plot, (b) The time variation of the action l x of the distribution of 
orbits. 

Figure 3: The numerical estimated local diffusion coefficients for 1\ = 3.0 
and 7.0, corresponding to m = 4 and 6 respectively, for 9 = 60° and v = 1.8 

Figures 4a and b: The diffusion coefficient D(Ii), as a function of the 
perpendicular energy of the charged particle, I\ for different angles of prop- 
agation, with a = 0.2 and v — 1.8. The □ are analytical values and the A 
are numerical ones: (a) For 9 = 20° and 35°. (b) For 9 = 45° and 60°. 
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